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Abstract. 

A unified electrodynamic approach to the guided-wave excitation theory is 
generaUzed to the waveguiding structures containing a hypothetical space-dispersive 
medium with drifting charge carriers possessing simultaneously elastic, piezoelectric 
and magnetic properties. Substantial features of our electrodynamic approach are: 
(i) the allowance for medium losses and (ii) the separation of potential fields peculiar 
to the slow quasi-static waves which propagate in such active media independently 
of the fast electromagnetic waves of curl nature. It is shown that the orthogonal 
complementary fields appearing inside the external source region are just associated 
with a contribution of the potential fields inherent in exciting sources. Taking account 
of medium losses converts the usual orthogonality relation into a novel form called the 
quasi-orthogonality relation. Development of the mode quasi-orthogonality relation 
and the equations of mode excitation is based on the generalized reciprocity relation 
(the extended Lorentz lemma) specially proved for this purpose with allowing for 
specific properties of the space-dispersive active media and separating the potential 
fields. The excitation equations turn out to be the same in form whatever waveguide 
filling, including both the time-dispersive (bianisotropic) and space-dispersive media. 
Specific properties of such media are refiected in a particular form of the normalizing 
coefficients for waveguide eigenmodes. It is found that the separation of potential 
fields reveals the fine structure of interaction between the exciting sources and mode 
eigenfields: in addition to the exciting currents (bulk and surface) interacting with 
the curl fields, the exciting charges (bulk and surface) and the double charge (surface 
dipole) layers appear to interact with the quasi-static potentials and the displacement 
currents, respectively. 
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1. Introduction 

Modern progress of material science and technology opens new potential possibilities 
in synthesizing complex and composite media with unique electromagnetic properties 
at microwaves and in optics. This requires revising some problems of guidcd-wave 
electrodynamics, in particular, the theory of waveguide excitation by external sources. 
The guided-wave excitation theory for passive media with isotropic, anisotropic and bi- 
anisotropic properties was developed in [1] . Elcctrodynamic processes in such passive 
media are fully described by Maxwell's equations and local constitutive relations, the 
most general form of which is inherent in bianisotropic media (BAM) with frequency- 
dependent constitutive parameters. Although magneto-electric phenomena in such 
media, by their microscopic nature, are brought about by non-locality of short- 
range polarization responses on electromagnetic actions [2-4], their macroscopic 
manifestations arc actually similar to those for real time-dispersive media. Indeed, 
the short-range character of these phenomena enables one to use the plane wave 
representation with wavenumber k = u/c so that in the first-order approximation 
all the constitutive tensor parameters of such a medium become solely local and 
frequency-dependent (see Ref. [2]). 

Unlike passive media, a true active medium requires, in addition to Maxwell's 
equations, an appropriate equation of motion for its elcctrodynamic description. 
Among active media we restrict our consideration to three kinds: 

(i) piezo-dielectrics with elastic properties providing the technological basis for 
acoustic- wave electronics [5 7] , 

(ii) dielectrics with ferrimagnetic properties (magnetized ferrites) forming the 
technological basis for spin- wave electronics [8 - 10] , 

(iii) nondegenerate plasmas with drifting charge carriers (in particular, semiconductors 
with negative differential mobility of hot electrons) constituting the technological basis 
for plasma- wave electronics [11-14]. 

For generality, we shall investigate a hypothetical spacc-dispcrsivc medium 
possessing simultaneously elastic, piezoelectric, ferrimagnetic and nondegenerate 
plasma properties in order to provide specific relations for any kind of the complex 
composite medium as a special case of the general situation developed below. 

Space-dispersive properties of such a waveguiding medium are related to non-local 
effects caused by specific interactions between adjacent particles of the active media, 
such as elastic interactions in piczo-dielcctrics, exchange interactions in ferrites and 
carrier diffusion effects in nondegenerate plasmas. Neglecting the non-local effects 
enables the equations of medium motion to be converted into the constitutive relations 
with frequency-dependent parameters. In other words, such a medium possesses 
properties of the usual time-dispersive active medium (TDAM) and in this case its 
elcctrodynamic properties fully conform to the ordinary anisotropic media examined 
in [1]. Our subsequent analysis will be based on results obtained there, among them a 
novel notion of the mode quasi-orthogonality for lossy waveguides (see section 3 of [1]). 

The objective in writing the paper is to generalize the guided-wave excitation 
theory developed previously for the waveguiding structures with time-dispersive 
bianisotropic media [1] to waveguides filled with a generalized (hypothetical) space- 
dispersive medium which contains all the piezoelectric, ferrimagnetic and plasma 
phenomena. Section 2 includes input information about the modal field expansions 
with separating potential fields to give a physical insight into the nature of the 
orthogonal complementary fields obtained previously [1]. Since the electrodynamic 
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description of space-dispersive phenomena requires, in addition to Maxwell's equation, 
special equations of medium motion, section 3 is devoted to consideration of the 
appropriate equations for three kinds of SDAM. Section 4 is fundamental and begins 
with an examination of the generalized reciprocity relation specially derived in the 
Appendix. This relation serves as a basis for obtaining the mode quasi-orthogonality 
relation and the equations of mode excitation taking into account specific contributions 
from the generalized hypothetical medium. Section 5 contains the general discussion of 
physical features in the mathematical description concerning the excitation of lossless 
and lossy systems valid for both BAM (see [1]) and SDAM waveguides. Mathematical 
notation here is the same as applied in [1], in particular, A stands for scalars, A for 
vectors, A for dyadics and A for tensors of rank more than two. 

2. Modal field expansions with separating potential fields 

The fundamental result of the previous examination [1] , which was obtained for passive 
media (including BAM as the most general case), is the incompleteness of cigcnmode 
basis for any waveguiding structure inside the excitation source region. This manifests 
itself in the fact that the modal field expansions Ea = X^fe ^kEk and Ha = '^f. AkH^ 
must be supplemented with the orthogonal complementary fields Ei, and which 
are related to the longitudinal components of the external currents and J™, i.e. 
the complete electromagnetic fields inside the source region are represented in the 
following form (cf equations (4.8) and (4.9) in [1]) 

E{rt,z)=Ea{rt,z) + Eb{rt,z) = J2Ak{z)Ek{rt,z) + Eb{rt,z) (2.1) 

k 

H{ru z) = Hain, z) + Jfb(rt, z) = ^ Ak{z) Hkin, z) + fffe(rt, z) (2.2) 

k 

where an unknown longitudinal dependence of the excitation amplitudes Ak{z) is due 
to exciting external sources as well as the complementary fields Ei, and H^. 

In [1] the complementary fields received a mathematical substantiation as the 
orthogonal complement to Hilbert space spanned by the eigenfield basis {Ek, Hk}, 
but their physical nature is still not properly understood. The present examination 
of space-dispersive media enables us to furnish an explanation of these fields as a part 
of the potential fields generated by external sources. 

The basic elcctrodynamic property of active media is associated with their ability 
to support propagating the special kind of slow waves (whose propagation velocity is 
much less than velocity of light characteristic of a medium under consideration) such as 
the surface acoustic waves (SAW) in elastic piezo-diclcctrics [5 7], the magnetostatic 
spin waves (MSW) in magnetized ferrites [8-10] and the space-charge waves (SCW) 
in semiconductors with negative differential mobility of electrons [12-16]. These 
slow waves, being of multimodal character for composite (multilayered) structures, 
constitute the so-called quasi-static part of the total electromagnetic spectrum, whose 
principal feature is related to the predominance of a relevant potential field (electric 
for SAW and SCW or magnetic for MSW) over its curl counterpart. 

Representation of the total fields {E, H) as a sum of their curl (£Jcj -f^c) and 
potential [Ep = — Vi^, Hp = — V^') parts is realizable on the basis of Helmholtz's 
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decomposition theorem [17]. In this theorem, the fields of every fcth eigenmode can 
be represented as 

Ek = Eck + Epk = Eck - Vipk Hk - Hck + Hpk = H^k - V^k (2.3) 

where V • Eck — and V • Hck = 0. 

Therefore, instead of one set of the total eigenfields {Ek, Hk} forming the basis of 
Hilbert space, we now have two sets involving the curl eigenfields {Eck, Hck} and the 
quasi-static eigenpotentials {^Pk, ^k}- This extends dimensionality of Hilbert space 
and allows us to anticipate the possibility of expanding the complementary fields Eb 
and Hb in terms of the scalar potential basis {(fk, V'fc}- 

Let us apply the vector curl- field basis {Eck, Hck} to expand the desired curl fields 

Ec{rt,z)=Y,Mz)Eck{rt,z) Hdn, z) Ak{z) Hckin, z) (2.4) 

k k 

and the scalar potential basis {(/Sfc, ipk} to expand the desired quasi-static potentials 
(p{rt,z) ^^Ak{z)ipk{rt,z) ipirt, z) Ak{z) ^pkin, z) . (2.5) 

k k 

In this case the complete fields inside the source region can be written, allowing 
for relations (pT^), in the following form 

E = Ec -Vip^^Ak (Eck - Vipk) - 2o X! 

k k 

= Y^AkEk - zoY^^^k (2.6) 

fc k 

H^ He - Vi^^Y.^^ (^-fc -'^^fc) - ^" II ^ 
fc fc 

= Y^AkHk- z^Y^'^i^k. (2.7) 



From comparison of equations (2.6) and (2.7) with (|2.1|) and (2.2) it follows that 



Eb = -ZoY^r'Pk Hb^-zoY^T^k- (2.8) 



az ^ — ' dz 

k k 



Formulae (2.8) give the expected expansions of complementary fields in terms of 
the quasi-static potentials of eigenmodes in a specific form involving the derivatives 
dAk/dz in place of the amplitudes Ak as expansion coefficients, which vanish outside 
the source region where Ak(z) = const. As evidently follows from (2.6) and (2.7), the 
complementary fields (|]^) are in fact a part of the total potential fields associated 
with external sources, which was formerly unexpandable, whereas the other part is 
included in the modal expansions. 

All the above results lead us to the important conclusion that the complemented 
Hilbert space spanned by two sets of base functions, consisting of the curl eigenfields 
{Eck, Hck} and the quasi-static eigenpotentials {^k, V'fcji is closed with respect 
to any function corresponding to arbitrary external sources because the desired 
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representations for the curl fields (2.4) and for the quasi-static potentials ( p.SD do 
not contain any orthogonal complements. So, if we entirely exclude the potential 
fields -Bp and Hp from our analysis, by using instead their scalar potentials (/? and V", 
together with the curl fields Ec and the appropriate sets of the mode quantities 
{(fik, 4'k} and {Eck, Hck} will constitute a complete basis that produces the modal 
expansions (2.4) and (|]^) with no orthogonal complements. In consequence, the latter 
fact provides the disappearance of the effective surface currents J^^f and Jg™f created 
by the complementary fields and given by formula (4.32) in [1]. 

I n th e subsequent examination, we shall consider that the modal expansions (2.4) 
and (2.5) involve the transverse distributions of physical quantities for all the modes 
of SDAM waveguide, which are taken to be known from a preliminary solution to the 
appropriate source-free boundary-value problem. Hence, the basic task is to find the 
mode excitation amplitudes Ak{z) inside the external source region. 



3. Constitutive relations and equations of motion for SDAM 

As was mentioned above, we consider the generalized (hypothetical) medium possess- 
ing simultaneously piezoelectrically-elastic, ferrimagnetic and plasma properties and 
restrict our consideration to the macroscopic model of SDAM, i.e. the medium under 
examination is regarded as a continuum characterized by pertinent phenomenological 
parameters. Such is the case in the long-wavelength approximation, when an excita- 
tion wavelength is much greater than typical medium dimensions such as interatomic 
distances for solids or Debye length for plasmas. For hot electrons in non-degenerate 
semiconductor plasmas, this holds true for the hydrodynamic, quasi-hydrodynamic 
and local-field approximations [14]. 



3.1. Piezoelectrically-elastic properties of a medium 

The stressed state of an elastic medium is specified by two second-rank tensors: the 
stress tensor T and the strain tensors'. The components Sij are related to components 
Ui of the vector u of medium particle displacement by the known relation [5 - 7] 

The components of stress tensor enter into the dynamic equation for the clastic 
medium written in the ordinary form of Newton's equation [5 - 7] 

where Pm is the mass density and U = du/dt = it is the medium particle velocity. 



The left part of equation (3.2) involves the total derivative with respect to time 
dU /dt = dU/dt + U • 'VU and the right part expresses the dynamic force F ='V -T 
exerted on a unit mass element (still without allowing for dissipative effects). 

If the elastic dielectric medium possesses piezoelectricity, both the strain S and the 
electric field E evoke an appearance of the electric polarization P and the elastic stress 
T. In this case the constitutive relations can be written in one of the conventional 
forms (ignoring magnetostrictive effects) [5 ~ 7] 

Pk = e^j S^j + eo xtkE, or P ^ § : S + to ■ E (3.3) 
T^J = cf^ki Ski - CHjEk or f=^^:S-i-E. (3.4) 



6 



In view of the relation D ~ e^E + P (3.3) can be rewritten as 

Dk = efcy S^j + ef^E, or D = ^ :S + e'-E. (3.5) 

In (3.3) - (^.5|), the quantities and e'^ = + X^) a-re the second-rank 

susceptibility and permittivity tensors, e is the third-rank piezoelectric stress tensor 
and is the fourth-rank elastic stiffness tensor (for brevity sake, the superscripts 5* 
and E will be dropped below). 

Acoustic- wave propagation losses in solids are cau sed by two dissipative effects 
which can be introduced into the dynamic equation ( |3.2D phenomenologically by means 
of the following quantities: 

~ f r . . 

(i) the internal friction stress T associated with the existence of viscous 
properties of an elastic medium, which should be added to the stress tensor T (to 
yield the total stress tensor = T + T^' ) in the form analogous to that used for an 
isotropic medium [7] 

Tl^^Tl^.M^ or T''^f,:S (3.6) 

where fj is the viscosity tensor considered as phenomenologically given; 

(ii) the dynamic friction force F^'' exerted by imperfections of a crystal lattice on 
the motion of acoustic phonons, which should be added to the total dynamic force 
F^ ~ V • in the form of a relaxation term 

F^' = -Tr^^p^U, or F'' = -t-' ■ p^-,U (3.7) 

where is the inverse rel axa tion ti me t ensor regarg ed a s phen ome nologically given. 



Allowing for relations (3^) and (3/7), equations {3A) and ( [3.2| ) are rewritten in 
the following form 



dS^j _ I (dUi dUj 



dt 2 \ drj dri 



(3.8) 



Pm^ = V•T- + Ff^ (3.9) 

In the case of pure harmonic processes (with time dependence in the form of 
exp(iti;t)), the above equations can be linearized in small-signal quantities (marked 
by subscript 1 unlike their static values marked by subscript 0) so that pm = 
PmO + Pml I U = Ui, etc with jp^il < IPmol- 



3.2. Ferrimagnetic properties of a medium 

Macroscopic dynamics of a ferrimagnetic medium uniformly magnetized by an external 
static field to the saturation magnetization Mq is described by the equation of 
motion written for the total magnetization vector M in the following form [8 - 10] 

dM 

— = - 7Aio (M X iJeff ) + R (3.10) 
at 

where 7 = jej/mo is the gyromagnetic ratio for magnetism of spin nature. 

The effective magnetic field H^s takes into consideration all the torque-producing 
contributions caused, in addition to Hq, by: (a) the Maxwellian field H (satisfying 
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Maxwell's equations), (b) the crystal anisotropy field = — Afc • M (due to 
magnetocrystalline anisotropy of a ferrite material), (c) the demagnetizing field 
Hd — — A/d • M (due to shape anisotropy of a ferrite sample), (d) the exchange 
field Hex ~ AexV^TW (due to nonuniform exchange interaction of precessing spins), 
namely [8-10] 

H,ff = + H - Af ■ M + AexV^M . (3.11) 

Here the net anisotropy tensor Af —Afc+Md allowing for both the magnetocrystalline 
anysotropy of a medium and the demagnetization anisotropy of a ferrite sample is 
assumed to be known, as well as the exchange constant Acx- 

The relaxation term R taking account of magnetic losses in ferrites is written in 
different forms, among them more convenient for us is the Gilbert form [8-10] 



R = a 



M dM\ 



with the damping parameter a considered as phenomenologically given or found from 
the resonance line half- width measurements as a = AH/Hq [10]. 

The Maxwellian field H, by its sense, is always a signal quantity, i.e. H = Hi, 
unlike the total magnetization M which is represented by separating a small-signal 
magnetization Mi in the for m M (r, t) = Mo + Mi(r, t) where [Mi] << \Mo\. Then 
the effective magnetic field ( ^.11 ) takes the following form 



H,s = Ho + Hi - AT ■ Ml + AcxV^Mi (3.13) 
where Hq = Hq — M • Mq is the static field inside a ferrite sample. 

3.3. Drifting charge carriers in a medium 

For the hydrodynamic description of nondegenerate plasmas with drifting streams of 
mobile charge carriers, it is more suitable and even necessary to apply, instead of the 
widely-used Eulerian description, a less known polarization description [13, 14, 18, 19]. 

As known [11-14], in the hydrodynamic model of non-degenerate plasmas the 
drifting charge carriers (say, electrons) are represented as a charged fiuid fiow 
characterized by such macroscopic quantities as the mean electron density n (or 
the charge density p = en), the mean electron velocity v (or the current density 
J = pv) and the electron temperature T (or the electron pressure p = nksT). 
Microscopic processes of scattering and thermal chaotic motion (or diffusion) of 
carriers are described in this model by such phenomenological parameters as the 
momentum relaxation time (or the mean time of free path) r and the thermal velocity 
Vt = {kBT/mY^^ (or the difffuson constant V ~ v%t). When intercarrier (electron- 
electron) collisions are rather frequent, there is local thermal equilibrium inside the 
carrier ensemble with the electron temperature T exceeding a lattice temperature Tq 
for high electric fields. 

In approximation of the local thermal equilibrium, the hydrodynamic force 
equation has the following form [11-14] 

^ + {v- V)v = —{E + V X B) ^ — ^—^ (3.14) 

ot m nm r 

where m is the effective electronic mass different from the mass mg of a free electron. 
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Equation ( |3.14 ) holds true for the case of uniform stationary heating of electrons 
when the electron viscosity, heat flow and thermal perturbations in an electron 
ensemble are negligibly small, so that there is no contribution from the so-called 
thermoforce and Vp = mwJVn [14]. Such a situation takes place when Tq << Tm, 
where Te is the energy relaxation time determining the rate of electron temperature 
perturbations and Tm = e/a = em/e^nT is the Maxwellian relaxation time determining 
the time scale of ac changes in the electric field and charge distribution. This condition 
means that the temperature keeps pace with signal perturbations in the electric field 
providing a local relationship between T and E [14]. The latter allows the momentum 
relaxation time r to be considered as a function of the electric field magnitude 
which is given phenomenologically or found from measuring the field dependence of 
mobility ^{E) = {e/m)T{E) and diffusion constant = v'^t{E). 

Strictly speaking, for plasmas placed in a magnetic field B the electron heating is 
produced by an electric field known as the effective heating field [14] 



Eh 



E^ + {b- Ef 



1 



where the quantity b = takes into account an influence of magnetic fields on the 
heating effect so that the quantities t, ^ and V now depend on Eh- 

For a small-signal situation when E = Eq + Ei, B = Bq + Bi, Eh = Eho + Ehi 
and all ac values (marked by subscript 1) are assumed to be much smaller than their 
dc counterparts (marked by subscript 0), we have 

rfr 



T{Eh) = T{Eho) + [j^] = To + Ti 



with To = T{Eho) and 



dlnr 



Ehl 



To \d\nEh J Ei^gEho 
Here we have denoted [14] 



(ko - 1) 



Ehl 
Eho 



(«o - 1) 



Eh 



Ei+Vq X Bi 
E'o 



(3.15) 



(3.16) 



where bo=/Lte-Bo and the anisotropy coefficient Ko = /Xd//ic (< 1, with inequality being 
true for hot electrons so that ko < for negative differential mobility) is defined as a 
ratio of the differential (/Zd) and static (/ie) mobilities which are equal to [14] 



(14 




Eq - 


h {bo • Eo)bo 




(1 + KO &0 ) + 


(1 + 




f (l-Ko) 




,-EoY/El 



dn{E)E 



dE 



and 



E=E. 



6 6 

He = KEho) = — T{Eho) = — To . 

m m 



Thus, the force equation ( 3.14| ) describes dynamics of the field-charge 
perturbations in hot electron gases characterized by the field dependence of momentum 
relaxation time t{E) given phenomenologically, which is known as the local field 
approximation [14]. 

The polarization (P) description holds an intermediate position between the well- 
known Lagrangian (L) and Eulerian (E) descriptions. To analyze the small-signal 
processes in drifting carrier streams it is customary to consider two states of the 
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stream - unperturbed (without a signal) and perturbed (with a signal). For the two 
states in the L-description the motion of a particular group of charges (located inside 
a physically infinitesimal volume called the liquid particle) is described by a time 
dependence of two radius- vectors ro{t) and r{t) appropriate to the unperturbed and 
perturbed stream. The fundamental dynamic variable of the polarization description 
is defined as a difference in these radius-vectors for two positions of the same liquid 
particle caused by signal action [14, 18, 19]: 

ri(ro,t) - r{t) - ro{t) (3.17) 

which is considered as a function of the unperturbed radius-vector rg and called the 



electron displacement vector. Hence, after introducing the displacement vector (3.17) 
the "trajectory" description of electron motion (typical for the L- variables) is replaced 
by the "field" description (typical for the E'- variables). Now we deal with a field of the 
electron displacement T'i(ro, t) which is completely identical to the field of the lattice 
particle displacement u{ro,t) applied in elasticity theory (see section 3.1). 

According to the Eulerian and polarization descriptions, the total instantaneous 



velocity v{r, t) of a particular group of charges satisfying the dynamic equation (3.14) 
is represented as [14, 18, 19] 

v{r,t) = Vo(r) + Mi(r,i) = ■Wo(''o) + Vi{rQ,t) (3.18) 

where the static velocity va is taken at two positions (perturbed, r, for the E- 
description and unperturbed, ro, for the P-description) of the same group of charges. 
The small-signal Eulerian (iti) and polarization (vi) velocities are defined by relation 



(3.18) at different space points but at the same point between them there is the 
following relation: Vi = Ui + {ri • 'V)vq. 

The polarization velocity Vi adheres to the equation of motion in the P- variables 



which is obtained from equation (3.14) in the following form [14, 19] 

— - + [vo ■ V)vi = 
at 



e 
m 



(^Ei + (ri • V)£;o + vi X Bo + vo X Bi + vo X (ri • V)Bo) 



N , V7 V7 ^ ^1 , Vo Ti + {ri • V)ro 
poV(V • ri) + Vri • Vpo \ (3.19) 

Po ^ ^ To tq To 



where the relaxation times tq and ri are given by equation (3.15). This equation allows 
for spatial nonuniformity in dc quantities Eq, Bq, po a-nd tq caused, for instance, by 
nonuniform doping of a semiconductor. 

Between the polarization variables Vi and ri there is the following rela- 
tion [14, 18, 19] 

vi = ^ + {vo- V)ri (3.20) 

ot 

which proves to play a role of the continuity equation 



(3.21) 
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Eulerian densities of charge pi and current Ji = PiVq + poUi are expressed in the 
P-variables in terms of the electronic polarization vector pi — p^Vx by the foUowing 
relations [14, 18, 19]: 

Pi = -V-pi (3.22) 
Ji = ^ + V X (pi X -yo) . (3.23) 



It is obvi ous that expressions ( 3.22 ) and ( 3.23|) satisfy the continuity 



equation (3.21) identically. The charge and current densities introduced by these 
formulae are fully the same as those produced by the motion of an actual dielectric 
with the polarization vector pi moving at velocity [20] . It is this fact that has given 
the name to the polarization description of mobile charges in vacuo and plasmas. 

3.4- Electromagnetic properties of a medium 

Electromagnetic fields are described by the usual Maxwell equations written for small- 
signal quantities as 

(3.24) 

V -01= pi V • Si = 

where the densities of charge pi and current Ji are represented in P-variables 



by (3.22) and (3.23). The induction vectors are associated with the polarization and 



magnetization vectors by the known relations [20] 

Di = cqEi + Pi and Bi = + Mi) . (3.25) 

For the generalized medium under examination, possessing also piezoelectric and 
magnetic properties, its polarization Pi and magnetization Mi adhere, respectively, 
to the constitutive relation (3.3) and to the equation of motion ( ^.1C| ) written in a 
linearized form. 

The polarization description of plasmas provides a convenient way to make an 
artificial replacement of plasma by an equivalent magneto-dielectric medium without 
mobile charges. Indeed, it is known that a real electric dipole pi moving at velocity 
Vq is perceived, by a fixed observer, as two immovable dipoles: electric pi and 
magnetic mi = pi x vq [20]. Hence, the polarization description operating with 
the electronic polarization vector pi allows any charged medium with mobile carriers 
to be represented as an equivalent polarized (with pi) and magnetized (with mi) 
medium with no mobile charges. In this case Maxwell's equations (3.24) by inserting 



expressions (3.22) and (3.23) can be written in the following form 



ot at 

(3.26) 

V • P)P = V • Bi = 

where the fundamental vectors Ei and Bi remain unchanged while the vectors Hi 
and Di are replaced with their equivalent polarization counterparts and D\. The 
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latter are introduced so that expressions ( [3.25 ) would hold their form with the use of 
these new vectors, namely 

Z)P = eo-Ei + and Bi = /xo(Hf + Mf ) (3.27) 

where the total polarization and magnetization vectors defined as 

= Pi + pi and Mf = Ml + mi (3.28) 

take into account the contributions from both a crystal lattice (Pi and Mi) and an 
electron ensemble (pi and mi = pi x Vo). Thus, the new equivalent field vectors 
entering into Maxwell's equations (3.26) in the polarization description are equal to 

= Di + pi and Hi + Va x pi . (3.29) 

Therefore, using Maxwell's equations of the "dielectric" form (3.26) enables one to 
consider all media including plasmas, as pure dielectric and describe them by applying 



the equivalent field vectors ( 3.29 ), which, in addition to the lattice polarizations 
(electric Pi and magnetic Mi), allow for the electronic polarization pi due to mobile 
charge carriers. In this case, the components of all field vectors, as follows from 
equations (3.26), are continuous on the boundaries of drifting carrier streams. Such a 
feature makes it more preferable to use the equivalent "dielectric" form of Maxwell's 
equations (without Ji and pi) instead of the usual form (3.24) (with Ji and pi). The 
point is that the latter generate the equivalent surface charge Pg'^ = n-pi and surface 
current J^"^ = (n • pi)vo on the boundaries of drifting carrier streams [14,15,21], 
which ensures discontinuity in the appropriate field components. 

4. Orthogonality and quasi-orthogonality of modes and equations of 
mode excitation 

The generalized reciprocity relation (the extended Lorentz lemma) 

Up to this point, we have considered the source- free region of a waveguiding structure 
with SDAM whose electromagnetic properties are described by the usual Maxwell 
equations (3.24) or their equivalent "dielectric" form (3.26). The latter form, which 
allows for the electronic polarization of drifting charge carriers, along with the lattice 
polarization and magnetization, is more suitable for subsequent examinations. Inside 
the region of exciting bulk sources (J^i, J^, p'^i, p™i) its curl equations are written 
for pure harmonic processes in the form 

V xEi = -iojBi - J^l and V x Hf = iojOf + J^i (4.1) 

where the electric and magnetic sources obey the continuity equations 

iojpli + V • J{;i = and iojp^i + V ■ J^^ ^ . (4.2) 



As follows from (3.27) - ( 3.29| ), in the equivalent magneto-dielectric medium 



without mobile charges, characterized by the total (lattice and electronic) polarization 
Pi = Pi+pi and magnetization Mf = Mi +mi, the field- intensity vectors Ei and 
= Hi + Vq X pi produce the flux-density vectors and Bi : 

^ eoEi + Pf = eaEi+ Pi + pi (4.3) 



Bi - poiHf + Mf) = fxo (Hi + Ml) 



(4.4) 
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where the vectors Pi , Mi and pi reflect physical properties of the medium discussed 
in section 3. 

A basis for deriving the equations of mode excitation is the reciprocity relation in 
complex conjugate form (often called the Lorentz lemma) extended to the generalized 
medium under consideration with piezoelectrically-elastic, ferrimagnetic and plasma 
properties. 

To derive the conjugate reciprocity relation it is necessary, in addition to the 
given system of equations (marked with subscript 1), to consider another system 
(marked with subscript 2 having also small-signal meaning) whose dynamic equations 
are all taken in complex conjugate form (see (5.25) and (5.26) of [1]). A conventional 
procedure applied to these equations gives relati on (5.27) (see section 5.2.1 in [1]), 
which can be rewritten with the aid of equations ( [4.3| ) and ( [l.4| ) in the following form 

V ■ (^Eix H^* + e;x Hf^ = - iuj(^Pi ■ E*2 - p; ■ Ei'^ 

- iuj)io(Mi • H; - M; ■ Hi) ~ iLu(pi • E'^ - pi, • E'l) 



Jb2 • -^1 



Jm* 
b2 



Hf 



(4.5) 



where E'12 — -£'1,2 + vq x Bi 2 is the electric field measured relative to an observer 
moving with the nonrelativistic velocity Vq. 

The first three terms on the right of equation (4.5) are calculated in the Appendix 
by the help of the appropriate equations of medium motion for elastic piezo-dielectrics, 
magnetized ferrites a nd dr ifting charg e carri er s treams (see section 3). Substitution 
of equations ( A. 10 ), ( A. 14 ) and ( A. 22 ) into (4J) finally gives the desired reciprocity 
relation (cf. equation (5.28) in [1]) 



V ■ S12+ qi2 = rjj^ . 
Unlike equation (5.28) of [1], here the total power quantity 



EM 



:1PM 



(4.6) 



(4.7) 



in addition to the usual electromagnetic contribution (with the polarization 
modification Hi 2 -fff,2 [13,21,22]) 



El X Hi* 



E*2 X Hi 



(4.8) 



contains also the contribution from non-electromagnetic fields in polarized media 



'-'12 — 



Tf . U*2 



T^* . Ui 



Tin* 
■ ''2 



*^2 



"^1 



■J^*+V2°*-J1 



(4.9) 



Here we have introduced new quantities: 

(i) for ferrimagnetic media - a vector J™ of the effective magnetization current 
density and a tensor V^™ of the effective magnetic (exchange) potential equal to 



J™ = iujfioMi 



(4.10) 
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(ii) for plasma media - a vector Ji of the electronic polarization current density 
and a tensor = V^^ + Vi^I of the effective electronic potential involving the 
electrokinetic potential tensor and the thermal (diffusion) potential equal to 



J^^ic^pi V{^^-Vovl v^'^^-^p, (4.11) 
e e po 

where — vi — (e/2m) (ri x Bq) is the resulting small-signal velocity of electrons in 
the polarization description allowing for their rotation in the static magnetic field Bq 
with the Larmor angular velocity ujl = — (e/2m)Bo [21, 22]. 

The relaxation processes in the generalized medium are taken into account by the 
sum of three contributions: 

qi2 = 2 {uj^ S;:^:Si+ p^o U; ■ t'' • (7i) + 2 ly,,^lo (^^^ (Mi • M*) 



(4.12) 



where we have introduced (i) for ferrimagnetic media the magnetic relaxation 
frequency i^m = c^^m = ck^PqMq and (ii) for plasma media the quantity 



^1 



Ti+ri- Vro 



(ko - 1) 



Eo Ei+Vq X Bi 



with the last expression obtained from ( ^.15| ) by ignoring spatial static non- 
uniformities (so that Vtq = 0) and assuming that the vector Bq is longitudinal (so 
that bg • Eq = 6oi?o and = Eq, as a result of (p^)) [14]. 

The interaction of the electromagnetic fields with the bulk external currents is 
taken into account by the following term in the right-hand side of (4.6): 



.(b) - 



Tin* 
^b2 



Hf 



(4.13) 



Representation of the electromagnetic field vectors {Ei^2 , -f^i 2) ^.s the sum of their 
curl (-Eci,c2 , -ffci C2) ^'^d potential (— Viy9i^2 , — VV'1,2) parts, as in (^.3D, gives 



V • (iJi X iff + m X hp) = V • [(£;ci X ifP* + e:, X /fP 



M^^B2y + ^;iiu;Bi) 



-{{iiopl.M + ikopl^rm) - {{^^Pl\W2 + (i^Pb2)*V'i) (4.14) 



where we have used the vector identity V- (Ax V0) = V'(0VxA) and equations ( [4.l| ) 
and 

The use of (4.14) in the reciprocity relation (4.6) leaves expressions (4.9) and ( [4.12D 
for 5'™ and qi2 unchanged, but modifies >S'f2'^ and which are now equal to 



K2 X 
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((pi(iwr>P)* + (^;(iwDP)) + (Vi(iwB2)* + ^2*(i'^-Bi)) (4.15) 



+ [{icjpl^M + (iiopl.r^i) + ((ic.pS\)^2* + (i^PS2)>i) • (4.16) 



Formula (4.K) involves three contributions to the transferred power from the 
electromagnetic curl fields, quasi-electrostatic potential fields and quasi-magnetostatic 
potential fields. Formula (4.16) reflects the interactions of the curl fields with the 
external bulk currents and the potential fields with the external bulk charges (surface 
exciting sources will be considered later). 

The generalized reciprocity relation (or the extended Lorentz lemma) in the 
form (4.6) (with equations (4.7) - (4.9 ), ([4.12D and (4.13) for the total electromagnetic 
fields or with equations (|4.7|), (4.9), ( [4.1 2] ), ( ^l5| ) and ( |4.16| ) for the fields separated 
into curl and potential parts) provides a basis for deriving the mode quasi- 
orthogonality relations and the equations of mode excitation in SDAM waveguides. 
To this end, the differential form (4.6) of the reciprocity relation should be converted 
into an integral form. 



.2. Quasi- orthogonality and orthogonality relations for SDAM waveguides 



Let subscripts 1 and 2 in the reciprocity relation ( |4.6| ) correspond to two eigenmodes 
with numbers k and I for a waveguiding structure with SDAM in the absence of 
sources, i.e. r^^-* = 0. Then integrating equation (4.6) over the total cross section S of 
the waveguide and applying the integral relation (2.25) in [1] give 

dPkiiz) 



dz 



+ Qki{z) = 



(4.17) 



where the complex cross-power (or self-power for k — I) flow transferred jointly by the 
kVa and /th modes is 



Pki(z) = Pt 



EM 
kl 



pPM _ 1 



jEM* 



Zq dS 



i("' 



EM 
kl 



AlAi e- (Tfc+^')^ 



Nki AIM e~ (-^i^ +-^')^ = - Nki alai (4.18) 



and the complex cross-power (or self-power for k — I) loss dissipated jointly by the 
fcth and /th modes is 



Qki{z) 



- I qUn,z) dS=- Mm AlAi e- ("^-^ +^')^ = - M^i 4 
•J s 



(4.19) 



The quantities^; , 5'™ and q^i in these formulas are obtained from (4.8) or 



( 4.15 ), (4.9) and ( 4.12 ) by replacing subscripts 1 an d 2 with k and /, respectively. The 
quantities Pm and Q^; defined by (4.18) and ( |4.19| ) determine: 
(i) the total power fiow (cf (2.34) and (2.46) of [1]) 



(4.20) 



k l>k 
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(ii) the total power loss (cf. equations (2.35) and (2.47) of paper [1]) 

Q(^) = E E Q''-' (^) = i E ^^fc (^) I' + ^ Re ^ E ^^'^^ «^ (^)«' (^) (4-21) 

fe i k k l>k 

where the mode amplitude is related to the excitation amplitude Ak by formula 
o-k — exp(— 7fez) for every fcth mode specified by the propagation constant 
Ik = CKfe + i/3fc and the set of cross-section eigenfunctions (marked with hat over 
them) {Ek{rt),Hk{rt), etc} such that Ek{rt,z) ^ Ek{rt) exp{-'ykz), Hk{rt, z) = 
Hk{rt) exp(-7fc 2:), e tc. 

In (4.18) - ( 1.21 ), following formulae (2.36) and (2.37) of [1], we have introduced 
the normalizing {Nki) and dissipative [Mki) coefficients constructed of the cross section 
eigenfield vectors. Expression for the normalizing coefficients {N^, 



kl 



EM 
kl 



given by formula ( 4.55 ) or ( 4.57 ), while the dissipative coefficients are equal to 



Mki = 2 



•J s . 



— 1 (M. 



■M, 



dS. 



(4.22) 



Formula (4.17) is in fact the required relation of mode quasi-orthogonality which 
reflects an independent transmission and dissipation of power by any one of mode 
pairs {k, I) in a lossy system (see section 3.1 of [ 1]). The expression of Pm and Qki 
in terms of iV^; and M^i is written in (4.18) and (4.19), owing to formulae (2.41) and 
(2.42) in [1]. Substitution of (4.18) and ( [4.19[ ) in equation ( [4.17D yields the desired 
form of the general quasi-orthogonality relation 



ilk +ii)Pki = Qki 



{l*k+li)Nki = Mm 



(4.23) 



which turns into the usual orthogonality relation for lossless waveguiding structures 
as a special case wi th Q ki — M^i = 0. 

Since relations (4.23) are completely coincident with the similar relations (3.6) and 
(3.7) of [1], the reasoning given there in section 3, concerning the quasi-orthogonality 
of eigenmodes in lossy waveguides (when Mki ^ 0) and the orthogonality of active 
and reactive eigenmodes in lossless waveguides (when M^i = 0), remains true for the 
waveguiding structures with SDAM. 



4-3. Equations of mode excitation for SDAM waveguides 



To derive the equations of mode excitation, the fields in the reciprocity relation (4.6) 
marked by subscript 1 (which will be dropped for exciting sources) are assumed to 
be the desired fields excited by the bulk and surface sources, whereas those marked 
by subscript 2 are the known fields of the fcth mode without sources. Any one of 
the physical quantities # (e.g. the components of electromagnetic fields, potentials, 
polarizations, etc) is written for the fcth mode as 

$k{ruz)^Sk{rt)e-^k^ (4.24) 

while the same quantity in the source region is represented in the complete form 

<Pi (n, z) = <Pa{rt, z) + <l'b(n, z)=^ Ai{z) <Pi{rt, z) + <Pb{rt, z) 
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= ^Ai(z)<l'z(rt)e-^'^+ <Z'fc(rt,z) = ^az(0)^i(rt)+ <?6(rt,2) (4.25) 

involving the modal expansion and the orthogonal complement The latter is 
true for all physical quantities except the curl fields and quasi-static potentials for 



which, according to equations (2.4) and (2.5), Ect = H^f, = and ipi, — Tpi, — 



The separation of potential fields causes the surface boundary conditions (formerly 
written for the total fields in the form of relations (4.6) and (4.7) in [1]) to be 
reformulatted. Inside the source region, besides the external bulk currents J^'™ and 
charges pj^'™, there are the surface sources, which are located along the contour Ls 
and given in the form of: 

(a) the current sheet with the electric and magnetic surface current densities 
J° and Jg™ which results in discontinuity in the tangential components of the curl 
magnetic and electric fields rig x H^^ and Ug x Ed, respectively; 

(b) the charge sheet with the electric and magnetic surface charge densities 
and p™ which results in discontinuity in the normal component of the electric and 
magnetic inductions Ug • D\ and ris • Bi, respectively; 

(c) the dipole (double charge) sheet with the electric and magnetic surface dipole 
densities J7g and 77™ which results in discontinuity in the quasi-static electric and 
magnetic potentials ipi and ipi, respectively. 

The corresponding boundary conditions written in pairs for the electric and 
magnetic sources have the following form [14,20]: 

(a) for the current sheets 

n+ X E+ + n- X E', = - J" (4.26) 
n+ X + n; X H^- = J° (4.27) 

(b) for the charge sheets 

nt-Df++n;.Dr= Pt (4.28) 
n+. B+ + n-.B- = (4.29) 

(c) for the dipole sheets 

n+ ^+ + n- = - r^l (4.30) 

nti>t+^s^i ^ (4.31) 
Mo 

Here, all the quantities with superscripts ^ are values taken at points of the contour 
Ls lying on its different sides marked by the inward unit normal vector . 

For deriving the equations of mode excitation it is necessary to integrate the 



reciprocity relation (4.6) (with replacing subscripts 2 by k) over the cross section 
S of a waveguide by using relation (2.25) in [1], with the result depending on whether 
the potential fields are separated. 

Without separating the potential fields (when Ei ~ Ea + Ei, and Hf — + 
in accordance with equation (4.25)), the substitution of equation (|4.8|) (with changing 



subscripts 2^k) into the integral relation (2.25) of [1] yields 
IS dzjs 



V . Sf^dS = / X iJP* + Elx HI) . z^dS 
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nb • LEb X Hi* + El X Hl]dl - 



El X Hi* +Elx Hi 



dz Js 



- ( J° • El + . //r) + ( J:=,ef • El + f • i?r) 



where 



5r = £;ix/fr + £^.*xiff 



(4.32) 



(4.33) 



Here we have used: (i) the boundary conditions (4.6) and (4.7) of paper [1] with the 
real surface currents J| and given on the contour with the unit normal vectors 
n^, (ii) the effective surface currents J^^f = — Tib x and J"],^ = rib x Et defined 
on the boundary contour Lb of the bulk source area S'b with the outward unit normal 
vector rib (see (4.24) - (4.28) in [1]). 

With separating the potential fields (when Ecb = H^^ — and ipb — i'b — ^ 
accordance with equations ( |2.4| ) and ( ^.5[ )), the substitution of ( 4.15| ) (with changing 
subscripts 2^k) into the integral relation (2.25) of [1] yields 



d 



S Jg 



E,i X Hi* + E:^ X HI, 



((^i(ic.r>P)*+^*(ic.r>P)) + (7^i(ic.B,)*+V*(ic.Bi) 



zo dS 



n+. [E,i X HI*+E:^ X ffP j + n; ■ [E,i x Hl*+ E*,, x iJP 
n+- (iPiiiioBk)* + rki'^ivBi)y + n; ■ (V'i(iwBfc)*+ ^^(iwBi 



dl 



dl 



dl. 



The use of the boundary conditions (4.26) - (4.31) for rearranging line integrals 
in the last formula gives 



:.E:,+jr-Hi*]di 



where 



Mo 



dl 



(4.34) 



Sfi' = [e^i X HI* + E:, X HI,) + {^liiiuDir + ^liiuDi: 
+ (M^^Bk)* + ijliiLoBi) 



(4.35) 



18 



It should be noted that, as distinct from equation (4.32), expression (4.34) does 
not contain the effective surface sources hke J^'^ since the curl fields (2.4) and the 
quasi-static potentials (2.5) have no orthogonal complements. 

Substitution of equation (4.9) (with changing subscripts 2^ k) into the integral 
relation (2.25) of [1] yields 



d_ 

dz 



Js 



PM 
fc ■ 



zo dS 



(4.36) 



where 



PM 



Vi ■ Jk* 



Tf * • Ui 



yr • Jk* + K^* ■ Ji 



(4.37) 



Integrating the reciprocity relation (4.6) (with changing subscripts 2^k) over the 



along with (4.36) and (4.37) give the following result: 

Qik{z) = Rik{z) 



cross section S' of a waveguiding structure and using ( [4.32 ), (4.33) or (4.34), (4.35) 

(4.38) 



dPifc(z) 



dz 



where 



s 



(4.39) 



Qlk 



qik dS 



(4.40) 



R,, = R^^,' + Rt^ = f r^^^ dS+ f 



dl. 



(4.41) 



The loss term Qik has the universal form determined by formula ( 4.12] ) (with 
changing subscripts 2— >fc) for qik- In contrast, the expressions for Pik and Rik are 
obtained different depending on whether the potential fields are se parated. 

Without se parat ing the potential fields, the power term ( 4.39 ) involves Sf^ and 
■^1 ™ giv en by ( |4.33| ) and (4.37), respectively. The bulk excitation term determined 
by ( 4.13 ) (with J^2™ = changing subscripts 2^fc) is equal to 



R 



(b) ^ 
Ik — 



Sb 



(4.42) 



and the surface excitation term follows from equation (4.32) in the form 



R 



Ik 



• HI 



dl 



(4.43) 



With separating the potential fields, the quantity Sf^ and Sf^ determining the 
power term (4.39) are given by (4.35) and (4.37), respectively. The bulk excitation 
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term determined by equations ( 4.16| ) (with J'^^ — p^'™ — and changing subscripts 
2 — > fc) is equal to 



I / • e * I • m / * \ 

+ [^^Ph fk + i^Pb Vk ) 



dS 



(AAA) 



and the surface excitation term, according to equation (4.34), take the form 



pis) 

^Ik 



Lb L 



Po 



(iujBk) 



dl. 



(4.45) 



The quadratic (power) quantities Pik and Qik in the form of ( 4.39| ) and (4.40) are 
constructed of the hnear quantities: (i) of the form (4.24) and (ii) — $a + 
of the form (4.25) involving the modal expansion 4>a = X^i^z^z and the orthogonal 
complement 4>i, (the latter is absent for the curl fields and quasi-static potentials). On 
this basis, we can rewrite ( 4.39| ) and ( 4.40| ) as 



Plk = Pak + Pbk 



;ik 



Qak + Q 



bk 



Pa 



Qak 



nEM 
^ak 



pPM 
^ak 



EM 
fc ■ 



Qak dS" . 



(4.46) 



(4.47) 



Here the subscripts a and b correspond to using the modal expansions #a and 
the orthogonal complements to construct the appropriate quadratic quantities; 
therewith Pbk = -P™ + = P,™ since always P^^ = by virtue of the fact that 



without separating potential fields Sf^-ZQ — because of Eh — ZqE^ and — ZqH^, 



whereas with separating them Sf,^ = because of ( |2.4| ) and (|2 



The basic property of the orthogonal complement to be (quasi-)orthogonal in power 
sense with respect to the fields of every fcth eigenmode, assures the following quasi- 
orthogonality relation for lossy systems: 



dPbkjz) 
dz 



+ Qbk{z) = 



(4.48) 



which is written in agreement with the mode quasi-orthogonality relation (4.17) 



Hence, with allowing for (4.46) - (4.48) formula (4.38) takes the form 



dPak{z) 



+ Qak{z) = Rlk{z) 



(4.49) 



involving solely the modal expansions with no orthogonal complements. 

Substitution of the modal expansions into expressions (4.46) and (4.47) for Pak 
and Qak allows us to reveal their longitudinal dependence (without an explicit writing 
of transverse coordinates): 

Pak{z) = PT{z) + PTiz) ^ J {sTiz) + STiz)) ■ Zo dS 
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I Js I Js 

= I](^w^+^w^^)^i(-z)e~^''^*+^'^" EE (^Y,NkiAi{z)e-^''^ e-^>- (4.50) 



Qak{z) 



I qak{z)dS ^y^Ai{z) I qikdS^ 
Js I Js 

Y,Mz) f <z:.,d5^(5]A4;AKz)e-^'^) 

J ^ S 1 



-Ik-' 



(4.51) 



By inserting equations (4.50) and (4.51) in relation ( 4.49| ) and representing the 
exciting integrals Rik in the form 



?(b) 



Rlk - ( Rk 



(b) 



Rk e 



involving the wave factor exp(— 7^2) explicitly, we obtain 



dA, 



A, ye 



Rk^ R 



(b) 



R 



(4.52) 



The square bracket in equation (4.52) vanishes owing to the quasi-orthogonality 
relation (4.23) and the required equations of mode excitation take the following form: 
(i) for the excitation amplitudes Ai{z) 



E^-^e--=<(z) + i?r(z), 



fc = l,2,. 



(4.53) 



(ii) for the mode amplitudes ai{z) ~ Ai{z)e 



(4.54) 



These equations have the general structure applicable for both cases of separating 
and not separating the potential fields. The only distinction between them consists 
in different forms of the normalizing coefhcients Nki = ^1;'^ + N^i^ (electromagnetic 
and polarized-medium) and the exciting integrals R^ = R^^"^ + (bulk and surface) 
which depend on whether the potential fields are separated: 



(a) without separating the potential fields (when equations (4.42) and (4.43) are valid) 



eI xH^i +Ei xHl*]-zadS- 



Ui + Tf . U*k 



V( 



Zq dS 



(4.55) 



Rk - R'i'^ + Ri'^ = 



• Ek + Jr -HTjdS- [J^ . Ek + JT ■ hT) dl 



■E, 



J^cf Hk 



dl 



(4.56) 
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involving necessarily the effective surface currents and Jg" f which were developed 
in [1] (see formulae (4.24 - 4.32)). 

(b) with separating the potential fields (when equations (4.44) and (4.45) are valid) 



PM 



vr-Ji + vc-Jk 



ZQdS (4.57) 



dS 



- [j^ . El, + • HTk) + (i^P^ + icopT ru 



eo Mo 



dl 



(4.58) 



involving no effective surface currents, unlike the previous case. 

All the quantities Nki and involve the cross-section eigenfunctions (marked by 
the hat sign above them) and the z-dependence of the exciting integrals Rk{z) is due 
to that of the external bulk and surface sources. 

The separation of potential fields has revealed a fine structure of the interaction 
between the external sources and the mode eigenfields (curl and potential), which is 
demonstrated by relation ( 4.58 ): the ciurents (electric ^ and magnetic J™^) interact 

with the curl fields (electric E^cfc and magnetic J^^j.), whereas the charges (electric p'^ ^ 
and magnetic interact with the quasi-static potentials (electric ^p, and magnetic 

V'fc). Also, there is the interaction of the displacement currents (electric iw-O)^ and 
magnetic \u)Bk) with the double charge (dipole) layers (electric rjl and magnetic Jyf ), 
if any. Usually, the latter do not exist in real physical situations, but can be introduced 
by the equivalence principle as equivalent surface sources. 

It should be mentioned that the generalized theory of guided-wave interaction, 
allowing for the potental fields, much like the theory elaborated here, was first 
developed by the author [22]. Both results are in good agreement only for lossless 
waveguiding structures, although the problem of the orthogonal complements went 
unnoticed then. Allowance for losses was made on the basis of the bi-orthogonality 
relation (instead of the quasi-orthogonality relation, as it is done here) which was 
obtained by introducing a subsidiary boundary-value problem (called the associated 
problem) to change artificially the sign of a loss parameter. In so doing, the eigenmode 
norms for lossy waveguides have no power meaning and the results of mode excitation 
pose serious difficulties in physical interpretation. 
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5. General conclusions on the excitation theory for BAM and SDAM 
waveguides 



The final equations (4.53) or (4.54) of the waveguide excitation theory developed are 
completely identical, by their structure, with equations (5.47) and (5.48) obtained 
in [1] for the waveguides with bianisotropic media. They constitute an infinite set of 
coupled differential equations of the first order in the desired modal amplitudes Ak 
or Ofc excited by a given distribution of the external sources (bulk and surface) which 



enter into the exciting integrals i?^*^"* and R^'j^^ defined by formulae (4.56) and (4.58) 



This set of excitation equations can be rewritten in the matrix-operator form 

N-Z{z)^R(z) (5.1) 

where the matrix-operator N formed from the normalizing coefficients is hcrmitian 
{Nki = N^f.) and the column- vector R has Rk as components. The column- vector Z 
is composed from the elements Zk — da^/dz + ^kO'k = {dAk/dz) exp(— 7^2;). 

Generally, the coupled equations developed hold true for dissipative systems since 



the coupling coefRcient Nki defined by equation (4.55) or (4.57) determines, according 
to (4.18), the cross-power flow Pki for any one of mode pairs, which is the case for 
lossy waveguides. It is evident that these equations remain true for nondissipative 
systems as a special case. 



For lossless waveguiding structures, the general quasi-orthogonality relation (4.23) 
turns into equation (3.10) of paper [1] to produce there the orthonormalization 
relations (3.18) and (3.24) for the active and reactive modes, respectively. They can 
be written jointly in the combined form 

{NkSki - for active modes 

(5.2) 
NkSj^i - for reactive modes 

where the norms Nk possess such a property that for an active mode it is real- valued, 
whereas for a reactive mode it is complex- valued and Nk = (see section 3.2 in [1]). 



Substitution of equation (5.2) into (5.1) yields 



{Rk/Nk - for active modes 

(5.3) 
Rj^/Nj^ - for reactive modes 



As follows from the orthogonality relations (5.2), every eigenmode of a lossless 
waveguide is orthogonal to all the other modes, except for the only one in combination 
with fields of which it forms its own norm: the active fc-mode is non-orthogonal to 
itself but the reactive fc-mode is non-orthogonal to its own twin-conjugate fc-mode for 
which 7fe+7| = 0. The fields of such a mode, non-orthogonal to the fc-mode, enter into 
the exciting integrals and determine the power of interaction with external sources, 
so as to supply the given fc-mode independently of the others (see section 3.2 in [1]). 
This is true only if the exciting sources are fixed, which lies outside the context of 
self-consistent treatment usually applied in practice. 



The independent excitation of every mode by given sources described by (5.3) is 
inherent only in lossless waveguides. For a lossy waveguide this is not the case. Even 
for the fixed exciting sources there is a dissipative coupling among modes expressed 
by non-zero off-diagonal elements Nki of the normalizing coefficient matrix TV. Such a 
dissipative coupling implies that, unlike lossless waveguides, the given external sources 
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excite the total set of eigenmodes as a whole rather than every mode as a single. It is 
very important to realize that this mode coupling of dissipative character is apparent 
since it exists only inside sources. Outside them even the presence of nonzero off- 
diagonal coefficients Nki leaves all the eigenmodes uncoupled. Indeed, for the source- 
free region the right-hand side of equation ( [3.1| ) vanishes {Rk ~ 0) and, if Nki ^ 0, 
from here it follows that Zk = or Ak{z) = constant. 

Hence, the eigenmodes of a lossy waveguide, being linearly independent solutions to 
the appropriate boundary-value problem, remain uncoupled outside the source region. 
Under these conditions, any eigenmode has the constant value of amplitude that was 
gained from sources at the exit boundary of their existence region and propagates 
along the waveguide without coupling to other modes. However, the picture of power 
transfer is more involved. Every fc-mode, besides the self-power flow = Pkk, also 
carries the cross-power flows Pki in pairs, together with the other ^-modes, which were 
also excited inside the source region and outside have constant amplitudes. According 



to the quasi-orthogonality relation in the power form (4.17), any one of mode pairs 
has a certain value of the cross-power flow Pki, as well as the cross-power loss Qu- 
This fact is a physical manifestation of the power non-orthogonality (called the quasi- 
orthogonality) among eigenmodes outside the source region. It is the existence of Pki 
proportional to the normalizing coefficients Nki that results in the apparent coupling 
of dissipative character among modes inside the source region, which is described by 



equations ( 4.53| ) or ( 4.54 ) with the exciting sources on the right considered as fixed 



At first glance, it may seem that the infinite set of coupled equations like ( p.lD 
obtained for lossy systems reduces to certain mathematical difficulties which are 



absent for lossless systems described by the uncoupled equations like (5.3). However, 
the difference between them disappears in the self-consistent formulation of a wave 
problem. In this case, the external sources appearing in Rk themselves can be 
represented as series expansions in terms of eigenmodes of another waveguide, which 
makes these sources exciting for the waveguide under consideration. Then, the 



uncoupled equations for single modes, for example (5.3), turn into an infinite set 
of coupled equations (see Refs. [13,22,23] where this technique is described). Solving 
such coupled equations for lossless systems differs little in computational complexity 
from an analogous solution for lossy systems. In both cases, to obtain the total solution 
of coupled equations, allowing for interaction between all modes, is unrealizable in 
practice. Usually, the major effect of interaction is determined by coupling among 
finite number of modes. Strongly interacting modes can be separated by the coupled- 
mode technique [24], which enables an approximate solution to be obtained with a 
precision sufficient for practical applications. 

Thus, we have developed a unified treatment of the electrodynamic theory of 
guided- wave excitation by external sources, applied to any waveguide with composite 
and multilayered structures involving complex media, with both bi- anisotropic and 
space-dispersive properties. 

Appendix A. Derivation of the generalized reciprocity relation for 
space-dispersive active media 

Appendix A.l. Contribution from piezoelectrically- elastic properties of a medium 



Let us rewrite equations (3.4) - (3.9) for the first system with subscript 1 : 



Ti^ij = Cijki Si^ki — Skij Ei^k (A.l) 
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Di^k — Skij Si^ij + €ik Ei^i (A. 2) 

lUjpmUi^i = ' - pm[/l,i (A. 3) 

Ti^., = Ti., + Tf;,^- , Ti^- = (A.4) 

For the second system, the similar equations are obtained from equations (A.l) 



through (A. 5) by taking complex conjugation and replacing subscripts 1 with 2 (both 
have the small-signal meaning). The equations obtained in such a way are numbered 
as (A.l') through (A. 5'), but will not be written explicitly for the sake of brevity. 
In order to calculate the term 

iw(Pi • e; - p; • = i^(-Di -e;- d;-Ei 



in the right-hand side of equation ( |4.5| ) it is necessary to: multiply (A.l) and (A.l') by 
—iujS2ij and iuSi^ij, respectively; multiply (A. 2) and (A. 2') by —iujE^i^ and iujEi^k, 
respectively; and to add all the results. Then 

ic. (Pi . e; -p;.e,) = xlo (ti,, s;,^ - t^,^ (a.6) 

where the relations eik = Cfei and Cijki — Ckiij have been applied. 
Multiplying (A. 5) by Tj*^^- and (A. 5') by Tiaj gives the following 

•^2.1 



V • (ti ■ u; + t; • Ui) + (ui. 



drj drj 



Multiplying (A. 3) by U^^i and (A.3') by C/i,i gives the following 

= 2 AnT^'C/i,^ = 2 p„,C/2* • • C7i (A., 
where the last equality is obtained by using the relation t^^ = tJI^. 
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To calculate the last term in the right-hand side of (A. 7) it is necessary to multiply 
(A. 4) and (A. 4') by dU^i/drj and dUi,i/drj, respectively, and to add the results. 
Then the use of (A. 5) and (A. 5') yields 



ii^Vijki Si^ki 



dm 



s: 



2,kl 



dUi, 
dri 



dU. 



2,1 



dU. 



2j 



dri 



- S. 



2,kl 



dri 



= Lo^rjijki (ySi^ij 82^1^1 + S'a^jj SiM^ — 2Lu'^r]ijki Si^ij S'j^^; = 2lu^ 82- r) : Si 
where the relations rjijki = rjkiij and rjijki = r]jiki have been applied. 



(A.9) 



By collecting formulae (A. 6) - (A.9) we finally obtain 



iujiPi ■ E; - P; ■ El] = - V • { Tf • U. 



2uj^ S;:fi : Si + 2 p,Jj; ■ t-'-Ui . 



(A.IO) 



Appendix A. 2. Contribution from ferrimagnetic properties of a medium 
For ferrimagnetic magnetic media it is necessary to calculate the term 

iuJ^io(^Ml ■ H*2 - m; ■ Hi'j = iuj(^Bi ■ h; - b; ■ h-^ 

on the right-hand side of (4.5). By using (3.1C) - ( 3.13| ) we write the linearized 



equation of motion for the first system marked by subscript 1: 



iujMi ~ — 7/10 

+ Aex Mo X V^Mi 



Ml X Ha + Mo X Hi - Mo X { Af ■ Mi) 
uj f Mo 



Mo 



X Ml 



(A.ll) 



where i>m = oluj^ — a7Mo-^^o is the magnetic relaxation frequency. 

The similar equation for the second system marked by subscript 2 in place of 1 
is obtained from equation (A.ll) by taking complex conjugation and numbered as 
(A.ll') (without its explicit writing for brevity). 

Vector-multiplying equation (A.ll) by M2 and taking account of the equalities 
Mi_2 • Ho = Mi_2 • Mo ~ valid for small signals we obtain 



;fMi X M^) = 



7Mo 



Mo (M2* • Hi^ - Ho (M2* • Mi^ - Mo (M2* ■ Af-Mi 



+ X,^Mo[M; -V'Mi 
The similar equation obtained in a such way from equation (A.ll') is 



-^..,^^[M;.Ml 

LUf,j Mo 



(A.12) 



iuj ( M2* X Ml 



7Mo 



Mo (Mi • i?2*) - Ho (Mi • M2*) - Mo (Mi ■ Af • M^ 



+ AexMofMi- V'M*) - tiy,, — ^(Mi-M;) . (A.13) 

V / UJ^„ Mo V / 



ojm Mo 
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By adding these equations with using the vector-dyadic identity V • (VA • B) = 
B • V^A + (V X A) • (V X B) + VA ; V-B and allowing for the symmetry of tensor 
A/ we finally obtain 

iuj^io f Ml • H; - M* • Hi 



V • (Vi" 



A/ Mo 



\ 2 

■)0 



Ml -Ma) (A. 14) 



where the magnetization current Ji^2 = i'j-'MoMi^2 and the effective magnetic 
(exchange) potential Vi™ = — AcxViVfi.2 of a ferrimagnetic medium have been used. 

Appendix A. 3. Contribution of drifting charge carriers in a medium 

For plasmas with drifting charge carriers, it is necessary to calculate the term 

iw(pi • E!; - pI ■ E[^ = iuj[pi ■ [E* +vax B;) - pi ■ {Ei + Vq x Bi)) 



on the right-hand side of (4^). To this end, we apply ( 3.19 ) and ( 3.20 ) for the first 
system with subscript 1 written in the form 

iujvi + {vq • ^)vi = 

El + (ri • V)£;o + Vl X Bo + Vq X Bi + vo X (n • V)Bo 



e 
m 



Po 



PoV{V ■ ri) + Vri ■ Vpo 



■th ^ Vo Ti + (ri • V)ro 

To To To 



(A.15) 



Vl — iijjri + [vo ■ V)ri 



(A.16) 



and the similar equations for the second system with subscript 2 in place of 1 obtained 
from ( |A.15 ) and ( A.16| ) by taking complex conjugation and numbered as (A.15') and 
(A.16') (without expli cit writing for brevity). 

Multiplying ( A.15| ) and (A.15') by — iLi;(m/e)po^2 ^^^^ itj(m/e)po^i, respectively, 
and adding the results with a combination of terms give 



ioj{pi-E'* -pl.E[ 



lujpo — <^iuj[r2 • Vl + ri • V2 



r; • {vo • V)vi - n • {vo • V)v; )■ + iujpo r* • (n • V)Eo - ri • (r* • V)Eo 



icjpoj rl ■ {vi X Bo) - ri • {v^ x Bo) + • (vq x (n • V)i 
n • (vo X {r; ■ V)Bo) I + it^^T^I Pori ■ V(V • ri) - pon ■ V(V • r*) 



r; ■ (Vri • Vpo) - ri • {Vr^ ■ Vpo) 
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1 TO 
lUjpo 

To e 



ri-\v. 



^2 - vo I - r. 

To 



Ti + ri-Vrp 
To 



(A.17) 



Let us transform the right-hand side of (A.17), with the last terms aheady in the 
desired form. 

The first term is transformed by using V • Jo ~ 0, and the result of 
multiplying ( A. 16 ) and (A. 16') by — iu!(jn/e)pQV2 and iuj{m/e)poVi, respectively, 
into the following form 



^2 • {vo ■ V)t>i - ri • {vo ■ V)t>2 



iLupo — <^iuj[r2 • vi + ri ■ V2 



V • 



lUJ — Vo[pi •V2-P2-VI 

e V 



(A.18) 



The second term vanishes after using the following vector-dyadic identity 

{AxB)-{V xC) = B-{A- VC) A ■ {B ■ VC) (A.19) 
with A — ri, B = and C — Eq, since V x Ep = 0. 



The third term is rearranged by employing ( |A.16D , (A. 16'), V-i?o = 0, V-Jq = 0, 
and the identity 

Ax{B- VC) -B x{A- VC) = VC • (B x A) - (V • C){B x A) 
with A ~ ri, B — r2 and C — Bq. Therefore, 



iLopo 



r2 • {vi X Bo) - ri • (^2 x Bq) + • (■Uq x (ri • V)Bo 



r,.(vo X (r2*-V)i?o)]| 



Bo • r2 X Vi - ri X 



^2) 



V- 



ILOVq 



(pi • {r; X Bo)/2 ~ p*2 ■ (ri x Bo)/2 



Vo • [VBo ■ (rj X ri 



(A.20) 



The fourth term uses identity ( [A.19| ) with A = ri, B = r2, C — Vpo and the 
fact that V x Vpo — 0, to be transformed into the following form 



\UJV. 



Por*2 ■ V(V • ri) - pori ■ V(V • r*)J + ^r* ■ (Vri • Vpo, 
ri • (Vr2* • Vpo)] I = '^ujvl ^ V • (por*2 V • ri - po^-i V • r*) 
(V • pi V • r2* - V • V • ri j + (r* • Vri - n • Vr2*) • Vpo 



V- 



e po 



P2 V • pi pi V • 



V • ( J2" + V2'" -^r 



(A.21) 



where the thermal (diffusion) potential — (fcsT/e)(pi^2/po) has been used. 
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Substituting ([A.18|) - (|A.2lD into ([A.17|), we finally obtain 



iu;(p, . E'* -pl-E[] - V . [V^ ■ Jr + Vr-Jl 



1 



Ti + r-i • Vrp 

To 



t»2 - 



Vro 



To 



(A.22) 



where the electronic polarization current J12 ^ i^Pi,2 and the effective electronic 
potential = ^i*;2 + ^1,2 ^ — (m/e) (i'o'i'i,2 + {''^t/Po} Pi, 2 ^ ) have been introduced. 



References 



[1] Barybin A A 1998 Modal expansions and orthogonal complements in the theory of complex 
media waveguide excitation by external sources for isotropic, anisotropic and bianisotropic 
media Progress In Electromagnetics Research vol 19 ed J A Kong (Cambridge MA: EMW) 
pp 241-300 

[2] Landau L D and Lifshitz E M 1960 Electrodynamics of Continuous Media (Reading MA: 
Addison- Wesley) 

[3] Post E J 1962 Formal Structure of Electromagnetics (Amsterdam: North-Holland) 
[4] O'Dell T H 1970 The Electrodynamics of Magneto-Electric Media (Amsterdam: North-Holland) 
[5] Auld B A 1973 Acoustic Fields and Waves in Solids (New York: Wiley) 
[6] Dieulesaint E and Royer D 1980 Elastic Waves in Solids (New York: Wiley) 
[7] Kino G S 1987 Acoustic Waves: Devices, Imaging and Analog Signal Processing (Englewood 
Cliffs NJ: Prentice-Hall) 

[8] Sodha M S and Srivastava N C 1981 Microwave Propagation in Ferrimagnetics (New York: 
Plenum) 

[9] Soohoo R F 1985 Microwave Magnetics (New York: Harper and Row) 
[10] Standi D D 1993 Theory of Magnetostatic Waves (New York: Springer) 
[11] Stix T H 1962 The Theory of Plasma Waves (New York: McGraw-Hill) 

[12] Steele M C and Vural B 1969 Wave Interactions in Solid State Plasmas (New York: McGraw- 
Hill) 

[13] Barybin A A 1977 Electrodynamic concepts of wave interactions in thin-film semiconductor 
structures Pt I Advances in Electronics and Electron Physics vol 44 ed L Marton (New York: 
Academic) pp 99-139; Pt II ibid vol 45 pp 1-38 

[14] Barybin A A 1986 Waves in Thin-Film Semiconductor Structures with Hot Electrons (Moscow: 
Nauka) (in Russian) 

[15] Barybin A A 1975 Boundary conditions on carrier stream surfaces in nondegenerate 

semiconductor plasmas J. Appl. Phys. 46 1684-96 
[16] Barybin A A 1975 Quasistatic solution of normal mode problem in semiconductor films without 

magnetic field J. Appl. Phys. 46 1697-706 
[17] Korn G A and Korn T M 1961 Mathematical Handbook for Scientists and Engineers (New York: 

McGraw-Hill) 

[18] Bobroff D L 1959 Independent space variables for small-signal electron beam analyses IRE Trans. 

Electron Devices 6 68-79 
[19] Barybin A A 1978 Mathematical treatment of polarization description of nondegenerate 

semiconductor plasmas Int. J. Electron 44 481-97 
[20] Panofsky W K H and Phillips M 1962 Classical Electricity and Magnetism (Reading MA: 

Addison- Wesley) 

[21] Bobroff D L, Haus H A and Kluver J W 1962 On the small signal power theorem of electron 

beams J. Appl. Phys. 33 2932-42 
[22] Barybin A A 1975 On the generalized theory of normal mode excitation in electromagnetic and 

polarized medium waveguides by external sources J. Appl. Phys. 46 1707—20 
[23] Barybin A A and Stepanova M G 1991 Theory of the coupling of optical modes in thin-film 

integrated-optics and acoustooptic structures Zh. Tech. Fiz. 61 120-6 (Engl, transl. 1991 

Sov. Phys. -Tech. Phys. 36 1147-50) 
[24] Louisell W H 1960 Coupled Mode and Parametric Electronics (New York: Wiley) 



